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I. 



INTRODUCTION 



A high frequency broad-band transducer which has a broad 
radiation pattern in a half-space is required for acoustic 
tracking of underwater vehicles in situations where the vehicle 
carries the sound source and the receivers are fixed. [1] The 
underwater ranges at the Naval Undersea Warfare Engineering 
Station (NUWES) , Keyport, Washington, are of this type. During 
tests of production torpedoes and new torpedo designs, the tor- 
pedo is tracked in three dimensions as it moves through the 
water . 

The tracking is accomplished using an array of hydrophones 
mounted on the bottom of the acoustic range which receives 
acoustic pulses radiated periodically from a projector mounted 
on each vehicle on the range. Currently, an acoustic projector 
mounted flush with the bottom of the torpedo sends out a 75kHz 
coded pulse that is picked up by the hydrophones. A shore- 
based computer receives the information from the range hydro- 
phones and computes the location and speed of the torpedo 
during the run. 

For the NUWES Keyport ranges, the signal to noise ratio re- 
quirements and the expected maximum transmission loss between 
source and the receivers places a minimum requirement on 
source level which is quite high. In addition, the pitch and 
roll motions of the torpedo and the effect of surface reflec- 
tion of the sound has led to a specification of a broad beam 
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width. Specifically, current requirements are for a uniform 
source level within a conical volume about 150° wide centered 
on the vertical axis of the transducer. 

The projectors are required to operate acoustically un- 
loaded (in air) 95 percent of the time, since the unit is 
turned on before launch and continues running after recovery. 
Another requirement placed on the projector to be used in 
higher speed vehicles is that the radiating surface be flush 
with the body of the vehicle. 

These various requirements place extreme demands on the 
physical properties of materials used in the transducer and 
create difficult problems of building a transducer which also 
has reasonable life expectancy and reasonable cost,. 

Another approach to this problem is the use of a flexural 
wave transducer as a broad-beam, broad-band -width source which 
would be readily adaptable to flush-mounting, a necessary re- 
quirement for high speed torpedoes. The possibility of using 
the torpedo hull itself as the radiating surface is especially 
attractive for the low drag, laminar flow shapes expected to 
become used in the future. Present designs use a transducer 
which is inserted into a recess in the hull. Although this is 
relatively smooth, it is not adequately smooth for ultra-low 
drag laminar flow. 

Some early work in this area was done by Barlow at NUWES 
(private communication) . The first previous work at the Naval 
Postgraduate School on this subject was done by LT Sevdik [1] . 
His calculations were based on classical plate theory as de- 
scribed in Morse and Ingard [2] and Malechi [3] , which does 
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not take into account of effects of transverse shear and 
rotatory inertia. Among the results of Sevdik's work was the 
noting of damping of the plate natural frequencies by fluid 
loading and of inconsistencies in the change of flexural wave 
phase speed with frequency for transducers of differing 
thickness . 

Mindlin [4] includes these effects in a three-dimensional 
plate theory and shows that classical plate theory begins to 
differ markedly from three-dimensional theory at frequencies 
above that at which the flexural wavelength is about ten times 
the plate thickness. His theory gives flexural-wave speed 
values for a given ratio of plate thickness to flexural wave 
length which are lower than those from classical plate theory. 

Feit [5] considers the case of an infinite plate excited 
by a simple harmonic point source and which is fluid loaded on 
one side. He discusses the significant changes in the theoreti- 
cal radiation patterns that result at frequencies above and 
below the critical frequency, that frequency at which the 
flexural wave speed in the plate equals the speed of sound 
in the fluid. His work is based on the Mind 1 in-Timoshenko 
plate theory. 

An experimental investigation was conducted of the effects 
of fluid loading on axially symmetric flexural waves in a 
circular aluminum plate 25 in. in diameter and 5/16 in. thick 
by LT Jarvis [7] in 1977. Jarvis' results are in good agree- 
ment with the Mind 1 in-Timoshenko plate theory for the speed 
of the flexural wave as a function of frequency. Jarvis' re- 
sults clearly show a change in the radiation resistance as 
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frequency increases, as evidenced by the increased damping of 
the plate natural frequencies with increasing frequency. At 
frequencies below the critical frequency, the fluid adds iner- 
tia which causes a lower flexural wave speed. Above the critical 
frequency, the wave speed increases abruptly to about the value 
for the unloaded plate, indicating that the radiation reactance 
(inertial) becomes very much smaller above the critical fre- 
quency. 

Until the work of Alper and Magrab [6] , a "closed form" 
solution of the coupled vibration of a finite circular plate 
clamped in an infinite, rigid baffle had not been obtained for 
either classical plate theory or the Mindlin-Timeshenko plate 
theory . 

The first treatment of this problem was given by Lord Ray- 
leigh [8] Mindlin [10] andDeresciewicz [9] have shown that the 
normal modes and natural frequencies of a plate of finite 
thickness are different and more numerous in a given frequency 
band that are given by classical plate theory. 

Reisman and Greene [11] have studied the "in vacuo" response 
of clamped circular plates to suddenly applied loads using the 
Mindlin-Timoshenko plate theory. They found that except for 
the first three or four natural frequencies of very thin 
plates (diameter to thickness ratio 40 or more) there was no 
similarity between the natural frequency spectra of the Mindlin- 
Timoshenko theory and classical theory. The difference in 
transient responses were also significant. 

In the Alper-Magrab paper the acoustic pressure distribu- 
tion is obtained in a half space bounded by an infinite rigid 
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baffle containing a clamped circular plate which is driven 
from the "in vacuo" side by a simple harmonic force having an 
arbitrary spatial distribution. A closed form solution is 
obtained which is valid for both the near and far fields. The 
analysis uses oblate spheroidal coordinates to obtain the solu- 
tion of the wave equation in the fluid. These coordinates are 
chosen because a bounding surface of one of the coordinates is 
represented entirely by the circular plate, and a bounding sur- 
face of another coordinate is represented entirely by the in- 
finite baffle with a circular opening the same size as the plate 
diameter. The third coordinate is identical to one of the 
polar coordinates, the coordinates in which the plate motion 
is obtained. As a result, the otherwise fairly complicated 
coupled problem is reduced to two simpler boundary value pro- 
blems which can be solved separately and the solutions coupled 
by the continuity condition at their common boundary. The re- 
sulting solution for the acoustic field is given by an expansion 
in spheroidal wave functions. Pressure at any point in the 
field can be evaluated by inserting the value of the wave func- 
tions at the point. 

This is in contrast to previous known attempts to solve 
this problem, all of which have resulted in the need to perform 
difficult double (surface) integrations numerically with each 
integration valid for only one point in the field. Furthermore, 
all known previous work has been restricted to the classical 
plate theory. 

In the following pages, a brief review is presented of 
Alper ' s theoretical work as reported in his doctoral thesis [16] 
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and from the point of view of using the results as a design 
tool for a transducer. 

As part of the present thesis, experimental work was conducted 
as a test of the Alper-Magrab theory, using measurements of the 
radiated beam pattern from a transducer. These results are 
presented and compared with Alper's calculation. 

Finally, a modification of the Alper-Magrab computer pro- 
gram is included as an appendix in a form suitable for running 
on the Navy Fleet Numerical Weather Facility Computer (CDC- 
65 00- "HAL" ) . 
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II. ' THEORETICAL CONSIDERATIONS 



A. APPROACH 

An outline of the approach to the problem [16] of calculat- 
ing the sound field from the flexual vibration in a clamped 
edge plate is : 

(1) The wave equation governing the motion in the acoustic 
fluid is transformed into oblate spheroidal coordinates. 

(2) Using the Mindl in -Timoshenko plate theory, the eigen 
frequencies and eigen functions of the plate are determined 

for a clamped circular plate using polar cylindrical coordinates. 

(3) The solution to the forced vibration problem is ex- 
pressed as the collection of all possible combinations of the 
plate eigenfunctions and is introduced into the appropriate 
nonhomogeneous differential equation governing the plate motion 
which contains a term representing the surface force acting on 
the plate [the surface force is the sum of the applied force 
and the fluid pressure (in terms of velocity potential expressed 
by the means of spheroidal wave functions) ] . The resulting 
equations are then reduced to algebraic equations in the two 
sets of arbitrary coefficients (plate motion and fluid motion) 
by using the orthogonality property of the "in vacuo" eigen- 
functions of the plate. 

(4) Because the velocities of the plate and of the fluid 
must equal each other at the interface of the plate on fluid, 
the appropriate derivatives of the plate displacement and fluid 
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velocity potential are equated. Orthogonality of the spheroidal 
^ave functions then yields another set of algebraic equations. 



(5) Combining (3) and (4) and solving the result gives 
equations which express the pressure at every point in terms of 
spheroidal wave functions, thus leading to a knowledge of pres- 
sure at the point of interest at a given frequency. 

B. SOLUTION 

The following steps are from the work of Alper [16] . As 
seen in Figure 1 a circular plate is clamped in an infinite rigid 
baffle, represented by the x, y plane. The upper part is filled 
with fluid and there is vacuum in the lower part. A spatially 
arbitrary, harmonic, transverse force is applied to the "in- 
vacuo" side of the plate. The objective is to calculate the far 
field pressure in the acoustic fluid produced by the forced motion 
of the plate. The equation of motion in the fluid is 



a 2 V 2 =IL 2 S - 2 | 



( 1 ) 



where 



Zl = ka, frequency parameter 



T = u,t, dimensionless time. 



0 = dimensionless velocity potential 



p - SI 2 J>0 dimensionless pressure 



( 2 ) 






V = -O.V0 dimensionless velocity 



dimensionless velocity 



(3) 



We will look for harmonic solutions of the form 




(4) 



where 
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0 = magnitude of velocity potential and is a function of 
space coordinates only. 

Hence, putting Eg. (4) into Eq. (2) gives 

p = ifl? 0 e 1 ^ = pe 1 ^ (5) 

s\ ■ T. j’ 

Again, p = iJL£) is a function of space coordinates only. 
Also, putting Eq. (4) into Eq. (1) gives 

a 2 V 2 (0e iT ) - il 2 (i 2 $ e iZ ) = 0 

A 

(a 2 tf 2 + -fL 2 ) 0 = 0 (6) 

Eq. (6) is the dimensionless Helmholtz equation. Its solutions 
must satisfy the boundary conditions in the fluid. Except for 
the conditions at the plate, these are: 

(a) The fluid velocity normal to the rigid baffle must be 
zero . 

(b) Only outgoing waves are accepted. 

Oblate spheroidal coordinates, which are shown in Figure 2 
[Van Buren, King, Baier, and Hanish, Ref. 17], are the most 
suitable for these boundary conditions. 

When the confocal ellipses and hyperbolas having an inter- 
focal distance of plate diameter 2a are rotated about the minor 
axis (z axis) , flattened ellipsoids and hyperboloids of one 
sheet are obtained. Under the transformation to dimensionless 
coordinates y and ^ . 



(7) 



z = a 



a [(l-7j 2 ) (\ 2 +l)] 1/2 


cos 0 


a [{I-T} 2 ) ( | 2 +1)] 1/2 


sin0 


a tyll. 





The constant surface is an ellipsoid of ra/olution. The 
cross section perpendicular to the x or y axis is an ellipse, 
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while that perpendicular to the z axis is a circle. In the xy 
plane the eccentricity of the ellipse is -i- and the semimajor 
axis is ^ . Hence the bigger ^ is, the larger the ellipse, 
and the more nearly circular it becomes. The surface^ = 0 is a 
circular disk of radius a in the xy plane. The surface 
is a sphere of radius °o . The = constant surface is a hyper- 
boloid of revolution of one sheet with anasymptotic cone whose 
generating line passes through the origin and is inclined at the 
angle cosV Q to the z axis. The cross section perpendicular to 
the z axis is an ellipse, sections perpendicular to the x or y 
axis are hyperbolas. The surface 7J= 0 is the xy plane except 
for the circular disk *£= 0. The surface 1^1 = 1 is the z axis. 
The intersection of ^ and constant surfaces is a circle. The 
angle ©, measured from the x axis towards the y axis determines 
a point on this circle. 

Domains are: 






n 



Eq. (23) becomes, in oblate spheroidal coordinates, 

.2. t 

^ c n ~ x 

0 = 0 



+ sL+nVrV) 



( 8 ) 



The general solution of this by separation of variables gives 



0 ,5 

/mn mn / rnn 



( 9 ) 
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where 



and 




cxDSmO / 

3in. <v\ G- 



Si^nC - !^^) and R mn Oi-f2.> i^) represent angular and radial 

functions respectively. Putting the general solution Eq. (9) 
into the Helmholtz equation Eq. (8) leads to two ordinary 
differential equations. 



d_ 

*1 



*1 



(■-$4 s i-^) 


+ 


A +-aY- m \] 

rv\n L 






. 7 7. 1 

A +il 5 _ m 


- 




-a l ^ 






( 10 ) 



I>° (u: 



where is the separation constant and m = 0,1, 2,3 

After mathematical calculations 

s4*n= fi w ©, ce) tl2) 

(n-m) are even values 

where A = constant, 

mn 



oo 



S mn tl ’ l - i M ) = 




mn 



S=0 



d C-ifl.) p ™ (71) 

1 l 



is the oblate angular function of the first kind, and 






is the associated Legendre function of order m and degree (m+Cj) . 
From [ref. 16, pp. 17-23]. 

Cl) 



C4) 

mn 



R ~ if' - R rnn <- in ' if > ‘ 1 R 1 if' 
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In the development of Eq. Cl 2 ) some of the possible solutions 
are rejected. The coefficient of the second independent solu- 
tion of Eq (4) , which is (-ill, ^ ) , is chosen to be 

zero, because at tj = 1 this function goes infinity. The coef- 
ficient of the other independent solution combination of Eq. (11) , 
which is (-i 12, ij ) /is also chosen to be zero, because 

only outgoing waves are considered. & R^'(-iil, i'S ) can 

mn mn £ 

be calculated in power expansion given by Flammer [ref. 18], 

Also by expanding Eq. (3) in oblate spheroidal coordinates yields 



V- 





^ I+l zd 7* 1 


. a? l\ 





(13) 



where , e and e are the unit vectors normal to the 

*1 \ e 

rl and Q surfaces, respectively. The velocity component normal 
to the Y surface is 



^ V '’iM 1 ^ 



(14) 



The baffle plane corresponds to 0. The rigid baffle 
condition requires that 



Vv 



= 0 = 



^=0 



n 



7=0 



Using Eq. (12) for Eq. (15) yields 

0 ) 



4 S 



rv\n 









= o 

7=o 



(15) 



(16) 
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After using the power expansion for (-iJ2 / V ) , odd 

values of (n-m) simply could not satisfy Eq . (16). Therefore 

odd values of (n-m) are rejected and hence it is found that 
Eq. (12) exists only for even values of (n-m) . 



Summation of all acceptable solutions yields 

oo 0° 



A 

t = 





( 1 ) 



A (-ill / ^)R^ ) (-in / if ) ®L(9) . (17) 

mn mn i mn L m 



m=0 n=m,m+2,m+4 



In the case of axisymmetric motion, m = 0 and then the 
far field solution becomes 



A. 

* 



of 



le £ 



oo 




"/a (1) 

‘- 1 * A on S o„ 



n=0 ,2,4 



(18) 



C. MINDLIN-TIMOSHENKO PLATE ANALYSIS [REFS 6 and 16] : 

The dimensionless equations of motion in polar coordinates 
for simple harmonic excitation are: 



K 



u. 1 15-/ ^ IFvl T 



=0 



(19a) 






(t + 2^jfL ^e_Vl 

\\9 u 7 - ae a*-/ a u se / 



= 0 (19b) 



23 



(19c) 



v u ^+$+ 



2 . 



Ca SL 

k T 



UT 



0^ 

h 



-^(u. ; 9)4-p a (a ; 9) 



T 



and 



(L^ + Mkw. 1_^HV 

T «■ ^/"S9 



Su 



where 



x , Z i 

il=U 



\jJ"{uj$)= dimensionless transverse displacement 

X = i^L 

U.,0 u ^u. u" 1 ^©- 2 " 

Y - Va\nr 

k = empirical shear coefficient for plate. 



'N 






^u,&) 

and are the angles of rotation of the nor- 
mal to the plate's middle surface a result of bending in the u 
and © direction respectively. 



u = 



_r 

CL 



(^(u,©-) = spatially distributed harmonic external load 
Pq.( u , 9)= reaction of the fluid pressure to £j (^U.,9^) 

A ^ 

4>u-H>r 



For 



1 



= 0 (free vibration) solutions to Eq. (19) are 



Tu- ^ v 1 U 



A 

29 



24 



( 20 ) 






in 

ae- 



A A A 

v\T = \N~\ + v\T^ , 

A A 



If , ur^ and 
^ ^ x A 

( v U/ e + c<, ) VAT 

< + < ) nO- 



A 

H satisfy 



- 0 

= 0 



— 7 1 

< \e + *i > H 



A 



= 0 



(21) 



where 
( <*, 



:e i ,Oi \ 

(V K ) 



( <T , <r ■ ) = 

11 K 




<*N ( S? i^-a 1 - ixl 

Solutions to Eq. (21) under the condition that they will 
be finite at u = 0 are: 

= 2 &,«, T„ (*, u.) ® 



*1 



m- O 






V>J 1- I ^2-M ^ C c ^2. u -') ®(V\ 

(V\=0 

W = zEL (f*3 U ) ^<V1 J 



SI 









( 22 ) 



rv\—o 



and 



DO 




av ^ 



\ 






*~s~ 



( 23 ) 
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Since J (ix) = i m I (x) only Eq. (22) can be used without 
m m 

losing generality. 

Boundary conditions for a clamped plate are 

ur(t,e)= o 

4^0,9) = o C24) 



Combining EQ. (20) and (22) yields 




Applying the boundary conditions, Eqs . (.24) to (2 5) and the re- 
quirement that solutions will be non-trivial the characteristic 
equations are obtained 



(■*,) Jm (*0 -(,(*»,) 




(26) 



where m = 0,1,2, ... and derivatives are w.r.t. u. The dimension- 
less frequencies , m= 0,1,2, ...; j=l, 2, ..., for which 

Eq. (26) is satisfied are the eigenfrequencies to which there 
are corresponding values < S mj . <r lmj and • The 

corresponding eigenfunctions are 
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A 

or. 
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T (V .1 X ( u a'lj© (&) 

J (V» V IMJ / "" T N ^ 1WJ / J ^ V 1 

J<v\ \ l_rv\t / ^ 



(a) 



m 0 / X/ 2 ^ «««/ j 1, 2 ^ 




m 0 ^ X/ 2 1 «•• / ^ * X , 2 ^ •■« * 

In the case of axi symmetric motion, m = 0, Eq . (26) takes 

the form (eiiminating the derivatives) 




and the eigenfunctions 



( 28 ) 
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A 

UT 

j 






% - - 1 A J . * ( V> ^ 



' u jT 

As. 

% 



= 0 



(29) 



because o. Thus, there is no coupling between U/" 

* A 

or and in axisyrnmetric notion. 

There can be rotations U) f 0 that are axisyrnmetric, but 

b A 

they will be independent of UT and (U . 

The characteristic equation for independent axisyrnmetric 

A 

Ve motion is obtained from Eq. (28) by setting the factor 

J| ( °y = °* (3o) 

When the roots of Eq. (30) are inserted into ^ , the _0_ 

a J 

values are obtained for axisyrnmetric motion. The ap- 

lied force is normal to the plate here and these modes are of 
no interest. This can be seen easily after calculating a value 



for 



w 



= I +■ 



*3 






For a thick steel plate, for example, — = 0.2,Xl_'/> 37 and 
for a thin plate XL would be even greater. 

Finally, Eq. (28) for this problem takes the form 

( 1)*, J 0 ( « z ) J x (^,) - ( G I- X («; )J 1 (^ Z ) = 0 (31) 

By using Eq. (21a) the roots Jfl- of this equation can be ob- 
tained. Symmetric eig enfrequencies for — = 0.04 and — = 0.2 

cl 3 

are given in Table I of Ref. 16. 
h 

when ' / <_ 0 / 

28 



_a.v a k T 
x 



when 



(32) 



Then calling 

|*> 2 = -oC 2 when t ^ le characteristic Eq. (31) for 

axisymmetric motion becomes 

( 07 _\ ) o<, I o (p) Jj ( ) + (^_-l )£ To C C<, ) I 1 (f>) = o (33) 

and corresponding eigenfunctions are 




D. COUPLING OF ACOUSTIC FLUID AND PLATE 



The following developments are also taken from Alper and 
Magrab [6]. On the ^ = 0 surface, which represents the plate, 
boundary conditions are: The normal fluid velocity and the 
acoustic pressure on the fluid are equal to the transverse 
plate velocity and the pressure on the boundary of the plate, 
respectively. 

Using Eq. (30) , 



V, 






l= c 



V, 




J_ 

7 






(35) 



(36) 



The external force (applied force per unit area plus 
acoustic pressure at the plate boundary of fluid) is given as 




( 37 ) 
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Hence, if Eq. (37) is substituted into Eq. (19c), then, use of 



Eqs. (19) and (36) couples the motion of the plate and the 

£ A 

acoustic fluid, so that solutions for yJ and VJ~ can be found 

Th- 



under the specified force distribution L . When CJ = 0 and 

A A A. 

JL = it , , the eigenfunctions ( \jj~, , uj - . , UJ. • ) satisfying 

r*j ^3 T 

combined Eq. (19) and Eq . (37) could be found. 

Let the solution for the plate motion be expressed as expan- 
sions in its eigenfunctions : 

0\ ^0 J — I 0 J 
A 

> • W - 




'I 



(38) 



<v\-0 J=-\ 

Putting Eq. (33) into the linear and homogenous combination 
of Eq. (19) and (37) there is obtained: 



(z s ) 2 2 

tv\-0 T-l 4 ° 



= 0 



= 0 



v =.0 J 

2- OQ 



(4) Z Z UT M j 



CL A 

(39 > 



r*=0 J-l 

The orthogonality condition for a clamped circular Mindlin- 
Timoshenko plate is defined as 

A 1 A A. "1 

lLdud& 



/ZTT /t r .f A A * A •) A 



A A ~ 
iy\j r\\C 



= ^rv\n 0 ^"om) n 



(40) 
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where 




is the Kronecker delta and 




( 41 ) 



Using Eq. (27a) and performing appropriate integrations 
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in which 
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For the ax i symmetric case, m = 0. 
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/N 

Multiplying the first of Eq. (39) by ^ J ^ ^ ^ ^ / 

the second by ^ Ul . U. dud O’ , and the third by 
A 

uj^,j U. cLu. d 9 r adding the three equations and integrating over 
0^0^ 7_rr f 0 ^ and using the orthogonality condition 

Eq. (40) leads to 

v<ue t42 , 

0 o 



m = 0,1,2,...; j = 1,2 



A 



By expressing f(u,Q) in a Fourier series f(u,0) = 



9C 



2 



m-o 



F m (u) =■ 
m 



lit 



® W O)<J 0 - 



(43) 



0 + *omV 

and using Eq. (12) to express ^Z/ there results 
°° 2. ^ po t\ ) (*») 

= T W ® Ce)- i nj J J CnC C-'^) 1^°) ©» 



A 

1 Li'W 
<Y\=0 



ftt=uO 



(44) 



Combining Eq. (42), (43) and (44) and performing integration 

with respect to 0 yields 

Co *■ 



hi) b .u , C jrLZ - sl \) = 



Q^ 

k 



e « r ! j (§4 2 °)3 



n=M,ivt+2 J tvy 






( 45 ) 



m 0,1,2,...; t 1,2, ... . 
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where 



o) 



9m, nj - I SfJ-iv) 3 a a 



and 



'mj - J F^ (u) \X^j (u-) ^ • <iu. 



Expressions for g . are 
r rn,n] 
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For the axisymmetric case m = 0. 

Combining Eq. (12) and (36) and using the first of Eq. (38) 
for w yields 

40 \ I n ^0 

i 7 B . W . W =- J -7 ft s 

Z_l < v\T ^ Zj ^ ^ Vdl 1 V 



J-l n = 



(48 ) 

The functions ( - iJ?_,7j?) , (n-m) even, are orthogonal 

on the interval Q4^([-4l * T ^ e orthogonality condition is 
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( 49 ) 



(n-m) even; m = 0,1,2, ... . 



where 



T = 

1 Art yy 



mn 






.i 



0o 



(g +2 m) ! £ d q C-i Wj 



°l=<U/< 



(Z<j+1w\+0c|! 



CSO) 



(n-m) even. 

Multiplying Eq. (48) by , ( - ill, 9^) integrating over 
the interval 0 • s °l v i n 9 f° r A mn and using the definition 
Eq. (46) yields: 



Do 






m = 0,1,2, 



; n = m , m+ 2 , m + 4 . 



(51) 



Putting Eq. (51) into Eq. (45) 

'* ° w 

1 



± 4 ? ffi ^ s ^ 

cs )z ,v% 



tv\4-l ; rv\4-*g 



(52) 



\ ^ ^s) Hwt 



m 0, 1, 2 , ...; t “* 1, 2 , ... . 

This represents an infinite number of sets of infinitely 
long linear non-homogenous equations in the coefficients D m j . 

For each value of m, there is a set of equations, t = 1,2, ... . 
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The infinite set of equations, yields solutions for ; j = 1, 

2, ... which, when put Eq. (38) gives the displacements of the 
plate and, when put into Eq. (51) , gives the coefficients A mn * 

i_ T_ 

The coefficients A are then used to determine the m 1 " harmonic 

mn 

' +- V* 

of the potential field using Eq. (12) and thereby the rn har- 
monic of the pressure using Eq. (5) . 

Instead of using an infinite set of equations there is a 
necessity to use truncation with a resulting limitation on 
accuracy (N = number of terms in the truncated solutions) . Alper 
[16] found that for the far field pressure using N = 11 for thin 
(^ = 0.04) and N= 16 for thick (|r - n 

a a = 0.2) plates gives solutions 

which are accurate to within 0.2 percent. 

Alper calculated patterns for the first three resonant fre- 

h h 

quencies for thin (— = 0.04) and thick (— = 0.2) steel plates. 

o. a. 

These are presented in Figures 3 and 4. Calculated pressures 
in the far field on the plate axis as a function of dimension- 
less frequency for thin and thick plates are presented in Figures 
5 and 6. Physical properties for steel plates with water as the 

fluid medium are: k = 0.86, V = 0.3, 0.42, IVn = 0.132 

T C 5 f J 

and P 0 = 1 . 

As a result it was decided that the best choice for building 
a transducer to meet the desired specifications would be a thick 
plate. A sketch illustrating the form of plate displacement for 
the first (il = 1.9) and second (-Q- = 6.9) symmetric modes is 
presented in Figure 7. 

The computer programs used for the theoretical calculations, 
which were designed for a CEC 6400 computer, were obtained from 
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Dr. Alper. These have been modified for the CDC 6500 and are 
listed in Appendix A. A shortage of time prohibited actual use 
of these modified programs for design work. 
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III. EXPERIMENTAL PROCEDURES 



An experimental test of the calculations of Alper and 
Magrab was conducted. A small, clamped edge plate transducer 
was constructed with dimensions which corresponded to the pub- 
lished calculated results for the thick plate, as shown in Fig- 
ure 8. The material used is steel. The ratio of plate thick- 
ness to radius, h/a, was chosen to be 0.2. The thickness was 
chosen so that the second, in-fluid resonance mode would be about 
75kHz. This corresponds to a value of the parameter = 6.9. 
Exterior dimensions were determined by the positioning and 
mounting hole on the cylindrically-shaped torpedo extension sec- 
tion, which served as the acoustic baffle for the transducer. 

The height of the outer cylindrical section at the rear was 
chosen to be about 3/4 wave length in order to create a high 
acoustic impedance for the plate edge in the neighborhood of 
75kHz . 

The driver for this transducer consisted of a piezoelectric 
ceramic longitudinal vibrator, similar to the sketch shown in 
Figure 9. The ceramic, a lead-zirconate titanate type, con- 
sisted of two axially polarized circular discs with a hole in the 
center. Masses on the ends and a clamping screw were constructed 
and their dimensions were selected to make the half-wavelength 
(free-free) longitudinal resonance to occur near 75kHz. A 
silicone grease was used between joints to ensure good acoustic 
coupling. 
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In order to provide a more nearly point source of excitation, 
the area of contact of the driver with the plate was reduced 
by the use of a small nut on the clamping screw between the 
driver and the plate. Later, measurements were made with the 
full face of the vibrator clamped to the plate. 

Resonance frequencies of the completed transducer assembly 
were determined using measurements of the electrical admittance 
of the piezoelectric elements. A Dranetz Impedance-Admittance 
meter was used for this purpose. Typical graphical results of 
conductance versus frequency are shown in Figures 12 through 14. 

The construction of the longitudinal vibrator used in the 
acoustic measurements prohibited the obtaining of a satisfactory 
isolation of the electrodes from electrical ground of the system, 
needed for making electrical admittance measurements with the 
unit in water. Therefore, resonance frequencies were not mea- 
sured in water. However, measurements were conducted of reson- 
ances of the driver with the transducer in air and in water 
using a different, isolated driver. This driver was a radially 
polarized cylinder of barium titanate, 2.54 cm long by 1.27 cm 
outside diameter. It was clamped to the plate using a screw 
and an end mass. Although the free-free half wave length 
resonance frequency for this element differed by several kilo- 
hertz from that of the other driver, the observed resonances 
of the combined systems differed only slightly. The change in 
resonance frequency caused by the immersion in water was 
within the uncertainty of frequency setting. The effects of 
radiation damping on the sharpness resonance was noted. 
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The acoustical measurements consisted of measuring the far 
field acoustic pressure as a function of angle. The arrange- 
ment of apparatus is shown in schematic form in Figures 10-11. 

A pulse method was used. Although the walls and floor of the 
water tanks are essentially anechoic, the reflection from the 
water-air surface can be significant, necessitating the use of 
a box-car integrator. The gate of the box-car integrator was 
adjusted so that only the sound wave traveling by a direct 
path from the source to the hydrophone was measured. 

The cylindrical baffle consists of an aluminum alloy cast- 
ing about 26.25 cm long and 32.4 cm outside diameter. Plastic 
end caps were used. The wall thickness is not uniform but is 
about 2.2 cm. A sketch is shown in Figure 12. 

The transducer is bolted into a hole which is located as 
shown. Lead weights were added inside the cylinder to make 
the assembly negatively buoyant. 

The rotation system could be arranged to measure the radia- 
tion pattern about either the longitudinal axis of the cylinder 
or the perpendicular axis. Measurements about both axes were 
made. The mechanical drive used to rotate the transducer in- 
cludes a potentiometer on the gear train which produces a D.C. 
voltage proportional to angle of rotation. This voltage is 
applied to the x axis of x-y plotter. The calibrated LC-10 
hydrophone used as a receiver and the transducer were 
separated by a distance of 1 m with the hydrophone and the 
center of the transducer at a depth of 1.92 . The 
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output of the hydrophone was amplified, rectified by the envelope 
detector and fed to the box -car integrator. In order to elimi- 
nate the effects of surface reflections, and minimize the ef- 
fects of standing waves in the tank a pulsed signal was used. The 
box-car integrator, gated to the direct path pulse produces a 
D.C. voltage that is proportional to the average amplitude of 
the received pulse envelope. This D.C. voltage is fed to the 
logarithmic converter which generates the ordinate of the beam 
pattern. The attenuator is used to calibrate the y axis to 
provide a decibel scale. It is needed because of the non- 
linearity of the envelope detector. 

Each radiation pattern was calibrated by adjusting the x-y 
recorder and driving voltage to obtain at least half scale 
reading on the integrator, and then by the means of attenuator, 
decibel increments on the y scale were marked. 

Results of radiation patterns about the perpendicular axis 
and the cylindrical axis are shown in Figures 20 and 21 and 
Figures 22 and 23, respectively. These were made with a small 
area of coupling between the driver and the plate. Then, with 
the small brass nut removed from the coupling between driver 
and plate the area of contact was increased to that of the 
disk at the end of the driver. Rotation was about the vertical 
axis- Results are shown in Figures 24 through 26. 

The box-car integrator gate width was suitably maintained 
to ensure that the integration would be done at steady state 
region of received signal. 

After some experimental work was done it was noticed that 
bubbles produced in the water greatly affect the results and 
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reduce loading effects of the fluid. Hence to prevent this 
effect a long brush was used for both hydrophone and transducer. 
Before each pattern was recorded, the source and receiver were 
brushed off in order that no effects due to bubbles would occur. 

A second series of measurements using the composite driver 
with small coupling area was conducted. The objective was two- 
fold. The first was to identify the modes of vibration associ- 
ated with the several resonance frequencies noted previously 
and, particularly, to determine which were the first and second 
modes. The second objective was to recheck the results of earli- 
er measurements. 

For the driver alone (free-free) conductance plotted as a 
function of frequency is shown in Fig. 13. Two resonances are 
noted. It seemed clear from observing the sensitivity to clamping 
provided by grasping it with one's fingers that the uppermost 
resonance involved flexural plate resonance of the lower and 
mass. The lower resonance appeared to be the usual free-free 
mode of vibration. 

When this driver is mounted on the transducer, a number 
of resonances can be observed. A plot of driver conductance 
as a function of frequency is presented in Figures 14 and 15, 
for the transducer alone and for the case where the transducer 
is mounted in the cylindrical baffle. 

A small capacitive displacement probe was used to measure 
the relative displacement amplitude at various positions on 
the outer face of the transducer. The probe consisted essen- 
tially of an 3 mm. diameter flat electrode which could be 
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positioned very close to the plate (about 0.1 mm) . This elec- 
trode and the plate itself formed an electrical capacitor. The 
electrical circuit provided polarization potential for this 
capacitor and measured the amplitude of the alternating poten- 
tial caused by changing capacitance when the plate vibrated. 
Although the probe measures the average effect of change in 
capacitance over a distance of 6 or 8 mm, it was sufficiently 
sensitive to get a reasonable estimate of the relative amplitude 
of plate motion as a function of distance from the center of the 
disk. A schematic representation of the displacement measure- 
ment apparatus is shown in Figure 16. 

A plot of the relative amplitudes is presented in Figures 
17 and 18 for each of the transducer resonances noted previously. 
The motion for the 01 and 02 modes was re-confirmed after the 
transducer was clamped in the cylindrical baffle. 
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IV. DISCUSSION OF RESULTS 



A. EFFECTS OF AXIS OF ROTATION ON PATTERN 

Radiation patterns were obtained for the same transducer 
but with different axis of rotation, one parallel to the cylin- 
der axis and the other perpendicular to the cylinder axis. The 
results for two such measurements are presented in Figures 20- 
21 and 22-23. 

It is concluded that although some differences between 
experimental patterns exist, these do not obscure the basic 
features of the radiation pattern to be used for comparison 
between theory and experiment. For purposes of convenience, 
all remaining patterns were made by rotating the baffle about 
its perpendicular axis . 

B. EFFECTS OF VARYING AREA OF CONTACT BETWEEN DRIVER AND PLATE 

Some radiation pattern measurements were made in which the 

size of the contact area between the longitudinal vibrates was 
increased from about 5 mm diameter to the diameter of the driver 
end mass, 17 mm. Since the diameter of the driver end mass is 
comparable to the dimension of the first nodal circle for the 02 
mode of the plate, it should be expected that the details of 
the displacement of the plate at the higher frequencies would 
not be a good approximation to that assumed for the point drive 
case. The radiation patterns made at high frequencies using 
this larger area of coupling were very complicated and difficult 
to interpret and they are not included here. The pattern at 
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low frequencies made with the larger area of coupling were not 
significantly different from those using the small coupling 
area. These are presented in Figures 24 through 26. 

C. IDENTIFICATION OF MODES OF VIBRATION OF THE PLATE 

The existence of a number of resonances in the neighbor- 
hood of frequencies calculated for the 01 and 02 modes of the 
transducer plate appear to be due to the frequency dependent 
impedance of the various parts attached to the plate. The 
driver element is one such attachment. However, the boundary 
conditions at the plate edge are probably more important. The 
massive cylinder at the rear of the plate was intended to 
provide a high mechanical impedance at the plate edge to simu- 
late the clamped edge condition. An outer flange of thickness 
comparable to (but larger than) that of the plate itself was 
needed to secure the transducer in the baffle. There were 
several modes of vibration which involved significant vibration 
of the outer flange. Some of these modes were damped where 
the transducer was clamped into the baffle. It is concluded 
that the resonances at 31 and 63.9 kHz were associated with 
the clamped edge plate modes having mode number 01 and 02, 
respectively. 

D. RADIATION PATTERNS 

The comparison between the radiation patterns calculated 
by Alper [16] and those measured in these experiments must take 
into account several factors which can give rise to differences. 

First, the assumed clamped edge boundary condition is only 
approximated in the tests. 
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Second, the baffle used is cylindrical in shape and is 
therefore neither infinite nor plane, so the assumption of the 
theory of an infinite plane baffle is not well approximated 
even though the diameter of the cylinder is about 7.5 times 
the diameter of the plate. 

Third, the fact that the baffle is not perfectly rigid per- 
mits some coupling of flexural wave energy from the plate into 
the baffle and the possibility of additional radiation of sound 
into the water from the baffle itself. 

Fourth, the area of contact between the driver and the center 
of the plate was made as small as possible but it was not a 
point driving force, as assumed in the theory. 

The comparison between theoretical and experimental radia- 
tion patterns should best be made by inspecting Figures 3 and 4 
and 27 and 29. The theoretical patterns due to Alper (shown in 
Figure 19) have been redrawn in the rectangular log-linear plot 
format used for the experimental data- 

The experimental pattern which should correspond most closely 
to the 01 mode is Figure 27 at a frequency of about 31 kHz. 

Those which should correspond most closely to the 02 mode 
are shown in Figures 20-22-29 at a frequency of 63 kHz. 

For the 31 kHz pattern it is seen that there are some small 
oscillations of about 3 dB near the center of a single major 
lobe. The level of this lobe at the 90° angle is down about 
18 dB below the axial value. This contrasts with the theoretical 
90° value of about -2 dB. The reasons for this disagreement 
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are not known, but it is believed that an important contribu- 
tor is the finite size and curvature of the baffle. 

For the 63 kHz patterns, there are some features which 
correspond well to the theoretical pattern. The angular 
position of the first minimum is at about 23 degrees, very 
close to the theoretical value of 22.5 degrees. The level of 
the measured maximum of the size lobe is about the same as that 
of the major lobe. The theoretical curves indicate a side 
lobe maximum level of -2.5 dB located at 48° and a level at 
90° of -5.6 dB . The experimentally observed level at 90® is 
about -20 dB. Again, this lower level near 90° is believed to 
be due largely to the finite baffle dimension. The additional 
structure on the measured patterns is believed to be due to 
radiation of sound from the baffle. 

Illustrations of the complexity of radiation patterns at 
other higher frequencies are shown in Figures 21-23-30. The 
frequencies here correspond to resonances in which there are 
additional displacements of the transducer parts particularly 
the outer flange (see Figures 16-17) . Very likely the baffle 
itself was also excited into some flexural modes of vibration. 

It is concluded that the major feature of the radiation 
patterns measured agree with those of the Alper theory and 
that the differences are probably due to boundary conditions of 
the plate which did not meet those assumed in the theory. 

E. RECOMMENDATIONS 

Since the goal for this effort was to develop the Alper- 
Magrab theory as a design tool for the development of a flexural 
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disk transducer that might be suitable for use as a sound 
source on underwater test range tracking, the following comments 
and recommendations seem to be appropriate. 

The adaptation of Alper's program to fit the CDC 6500 Com- 
puter, HAL, is virtually complete. Shortage of time precluded 
completion of this part of the work. It is recommended that the 
work be finished and that it be used as intended. 

It appears that some sort of variation on the radiation 
pattern from the second radially symmetric plate mode does offer 
the potential for controlling the source level near the axis 
and yet maintain a relatively broad beam width for the secondary 
lobes. The patterns observed here show a deep but narrow null 
at about 22 degrees off the central axis. It may be possible 
that this null could reduced by proper distribution of excitation 
force on the plate. It also is possible that the effect of 
this narrow null on range tracking is acceptable. It is recom- 
mended that the experience of this work be employed to con- 
struct and test on the tracking range a flexural disk trans- 
ducer which operates in this mode. 

The smoothest radiation pattern is given by the plate 
operating in the first symmetric mode. The transducer used in 
these tests was not optimized for this mode. It is likely 
that a transducer designed to operate in this mode at 75 kHz 
would have such a small plate diameter that a problem would 
exist in fitting an adequate sized driver into the available 
space. If the program development tool becomes workable, cal- 
culations should be made to test the design at 75kHz to deter- 
mine feasibility. 
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A much more likely application for a flexural disk 
source operating in the first mode appears to be at frequencies 
significantly lower than 75 kHz, such as are used on other test 
ranges. 
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Fluid Side 




Vacuum Side 



Figure 1, Geometry of Problem 
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Figure 2. Oblate Spheroidal Coordinate System 
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Figure 3. Far Field Beam Patterns for Thin Plate 
(From Alper , Ref. 16) 
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Figure 4. Far Field Beam Patterns for Thick Plate 
(From Alper, Ref. 16) 
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Figure 5» Pressure in Far Field on Plate Axis as a Function of 
, Dimensionless Frequency. (From Alper, Ref. 16) 
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Dimensionless Frequency. (From Alper , Ref. 16) 




Figure 7 . Form of Plate Displacement 
(Schematic) 
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Figure 8. Transducer Dimensions 
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Figure 9, Driver Elements and Position on the 

Transducer 
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Figure 10. Radiation Pattern Measurement Apparatus (Transmitting) 
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Figure 11. Radiation Pattern Measurement Apparatus 

(Receiving) 
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Figure 12. Cylindrical Baffle 
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Figure 13. Conductance/Frequency Plot 
Driver Alone (with nut) 
in Air 
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Figure 14. Conductance/Frequency Plot. Transducer in Air. 
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Figure 15. Gonductance/Frequency Plot. Transducer mounted in Baffle in Air. 
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Figure 16. Relative Plate Displacement Measurement Apparatus. 
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Figure 17. Relative Plate Displacement of Transducer. 
(Unmounted, in Air) 
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RELATIVE PLATE DISPLACEMENT 




RELATIVE PLATE DISPLACEMENT 




Figure 18. Relative Plate Displacement for 01 and 02 
mode (transducer mounted in baffle in air) 
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Figure 19. Theoretical Radiation Patterns of 01 and 

02 Mode. 



68 




Figure 20. Radiation Pattern (small area, perpendicular 

axis) 
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Figure 21 . Radiation Pattern (small area, perpendicular 

axis) 
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F£EQ.u£MCY= 69.5 kHi. 




Figure 22. Radiation Pattern (small area, perpendicular axis) 
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Figure 24, Radiation Pattern (large area, perpendicular axis) 
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FREQUENCE 31-57" ktk 
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Figure 28. Radiation Pattern (small area, perpendicular axis) 
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FREQ.u£N<iy= 92.93 




Figure 30. Radiation Pattern (small area, perpendicular axis) 
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APPENDIX A 



COMPUTER PROGRAMS FOR FAR FIELD PRESSURE 
AND ASSOCIATED PARAMETERS 

Alper ' s original programs were not changed very much. 

For the main program, the outputs of the Program, Roots, and 
the Program, Caller, (for getting sum') are needed. After 
getting the output of Program, Roots, Program, Caller output 
is obtained by using these results: Finally, the main pro- 

gram will be run by using the obtained parameters and sums. 

For this reason some Formats were changed to fit inputs and 
outputs each other. The only problem which still was not solved 
is that of Subroutine Bessel. 
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